term=4; n=400;

data(term);

load data;

A=matrixA(term,n);

b=source(term,n,uex{1}(0),uex{1}(1),uex{2}(0),uex{2}(1),uex{3}(0),uex{3}(1),uex{4}(0),uex{4}(1));

u=A\b;

maxL2err(u,n,term);